An RF-Driven Josephson Bifurcation Amplifier for Quantum Measurements 
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We have constructed a new type of amplifier whose primary purpose is the readout of supercon- 
ducting quantum bits. It is based on the transition of an RF-driven Josephson junction between 
two distinct oscillation states near a dynamical bifurcation point. The main advantages of this new 
amplifier are speed, high-sensitivity, low back-action, and the absence of on-chip dissipation. Pulsed 
microwave reflection measurements on nanofabricated Al junctions show that actual devices attain 
the performance predicted by theory. 



Quantum measurements of solid-state systems, such as 
the readout of superconducting quantum bits 0, 0, 13, S 
0- E3- ■ challenge conventional low-noise amplification 
techniques. Ideally, the amplifier for a quantum mea- 
surement should minimally perturb the measured system 
while maintaining sufficient sensitivity to overcome the 
noise of subsequent elements in the amplification chain. 
Additionally, the characteristic drift of materials proper- 
ties in solid-state systems necessitates a fast acquisition 
rate to permit measurements in rapid succession. To 
meet these inherently conflicting requirements, we pro- 
pose to harness the sensitivity of a dynamical system - a 
single RF-driven Josephson tunnel junction - tuned near 
a bifurcation point. The superconducting tunnel junction 
is the only electronic dipolar circuit element whose non- 
linearity remains unchanged at arbitrary low tempera- 
tures. As the key component of present superconducting 
amplifiers [3, 111 LUj] , it is known to exhibit a high degree 
of stability. Moreover, all available degrees of freedom in 
the dynamical system participate in information transfer 
and none contribute to unnecessary dissipation resulting 
in excess noise. The operation of our Josephson bifur- 
cation amplifier is represented schematically in Fig. 1. 
The central element is a Josephson junction whose criti- 
cal current Iq is modulated by the input signal using an 
application-specific coupling scheme (input port), such as 
a SQUID loop or a SSET 0. The junction is driven 
with an sinusoidal signal inp sin(arf) fed from a trans- 
mission line through a directional coupler (drive port). 
In the underdamped regime, when the drive frequency u 
is detuned form the natural oscillation frequency uj p and 
when the drive current Ig < Irf < Ib *C Iq, the system 
has two possible oscillation states which differ in ampli- 
tude and phase [ill Il2| . Starting in the lower ampli- 
tude state, at the bifurcation point inp = Is the system 
becomes infinitely sensitive, in absence of thermal and 
quantum fluctuations, to variations in Iq. The energy 
stored in the oscillation can always be made larger than 
thermal fluctuations by increasing the scale of Io, thus 
preserving sensitivity at finite temperature. The reflected 
component of the drive signal, measured through another 
transmission line connected to the coupler (output port), 
is a convenient signature of the junction oscillation state 
which carries with it information about the input signal. 
This arrangement minimizes the back-action of the am- 
plifier since the only fluctuations felt at its input port 



arise from the load impedance of the follower amplifier, 
which is physically separated from the junction via a 
transmission line of arbitrary length and can therefore 
be thermalized efficiently to base temperature. In this 
Letter, we present an experiment that demonstrates the 
principle of bifurcation amplification, and we show that 
the sensitivity obtained is in good agreement with a well- 
established theory of dynamical transitions, thus making 
it relevant for stringent tests of superconducting qubits 
and gates. 

The dynamics of the junction are descibed by the time 
evolution of the junction gauge invariant phase differ- 
ence S (t) = J_ dt'2eV(t')/h where V is the voltage 
across the junction. In the presence of a microwave drive 
irf sin(wi), the oscillations in the junction phase can 
be parametrized using in-phase and quadrature phase 
components S(t) = S\\ sin(uit) + S± cos(wi) (higher har- 
monics of oscillation are negligible) . When the detuning 
a = (1 — oj/uip) and the quality factor Q — ui p RC sat- 
isfy aQ > v3/2, then two steady state solutions can 
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exist for S(t). Here u p = (2eIo/hC) is the junc- 
tion plasma frequency, C is the capacitance shunting the 
Josephson element and R = 50 SI is the characteristic 
impedance of the transmission line at the output port. 
The two oscillation states appear as two points in the 
(Su, 5±) plane and are denoted by vectors labelled and 
1 (see Fig. 2). The error-current 01 which describes 
the generalized force felt by the system is also plotted in 
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FIG. 1: Schematic diagram of the Josephson bifurcation am- 
plifier. A junction with critical current Io coupled to the input 
port is driven by an RF signal. In the vicinity of the dynam- 
ical bifurcation point, the phase of the resulting oscillation 
depends critically on the input signal and manifests itself in 
the reflected signal phase <j). 
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(<5||,<5j_) plane. Its value goes to zero at the attractors 
corresponding to states and 1 and also at a third ex- 
ternum which is the dynamical saddle point. Also shown 
in Fig. 2 are the calculated escape trajectory from 
state (dashed) and the corresponding retrappping tra- 
jectory [Tj| into state 1 (solid line). Fig. 2 has been 
calculated for a = 0.122, Q = 20 and i RF /lB = 0.87, 
where I B = 16/(3^3) a 3 / 2 (l - af/ 2 h +0{l/{aQ) 2 ). 
These values correspond to typical operating conditions 
in our experiment. The dynamical switching from state 
to 1 is characterized by a phase shift given here by 

= — 139 deg. Using the 



200 



tan 



(*J-*fi)/(*i 



si 



junction phase- voltage relationship and the transmission 
line equations, we can calculate the steady-state magni- 
tude and phase of the reflected microwave drive signal. 
The change in the oscillation of S results in a shift of the 
reflected signal phase A<^>oi = 89 deg. Since there is no 
source of dissipation in the junction chip, there should be 
no change in the magnitude of the reflected signal power, 

even though y^jj - 5°) 2 + (S]_ - 5° ± ) 2 ^ 0. 

Our main sample consisted of a single shadow-mask 
evaporated AI/AI2O3/AI tunnel junction with Iq = 
1.17 (iA, shunted with an on-chip lithographic capaci- 
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FIG. 2: Poincare section of an RF-driven Josephson junction 
in the bistable regime [a, = (1 — oj/uip) = 0.122, irf/Ib = 
0.87). The coordinates <5m and 8± are the in-phase and 
quadrature-phase components of the junction gauge- invariant 
phase difference 8. The color code gives the magnitude of the 
error current i e which describes the "force" on S. The two 
stable oscillation states, labeled by and 1, are indicated by 
white line segments. The basins of attraction corresponding 
to the two states are separated by the blue dotted line (sep- 
aratrix). Point S which lies on the separatrix is the saddle 
point at which the escape trajectory from state (dashed 
line) meets the retrapping trajectory into state 1 (solid line). 
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FIG. 3: Reflected signal phase <j> as a function of drive cur- 
rent Irf /Iq- Symbols denote the mode of <f> with up and down 
triangles corresponding to increasing and decreasing current, 
respectively. The dotted line is (<fi). The calculated bifurca- 
tion points, Ig and Ib, are marked on the horizontal axis. The 
and 1 phase states are reminiscent of the superconducting 
and dissipative states of the DC current biased junction. 



tance C = 27.3 pF [3 to obtain a reduced plasma fre- 
quency bj p /2ii = 1.80 GHz. The dynamics of the transi- 
tion between the two oscillation states were probed using 
microwave pulses, generated by the amplitude modula- 
tion of a CW source with a phase-locked arbitrary wave- 
form generator with 1 ns resolution. The reflected sig- 
nal was passed through an isolator at base temperature 
T = 0.28 K to a matched HEMT amplifier at T = 4.2 K 
in the spirit of Day et al. 16]. At room temperature, 
the reflected signal was further amplified, mixed down to 
100 MHz and finally digitally demodulated using a 2 GS/s 
digitizer to extract the signal phase <j>. 

We first probed the drive current dependence of the 
reflected signal phase ^{irf) by applying a 4 /is long 
symmetric triangular shaped pulse with a peak value 
0.185 Jo- The demodulated reflected signal was divided 
into 20 ns sections, each yielding one measurement of (f> 
for a corresponding value of irf- The measurement was 
repeated 6 x 10 5 times to obtain a distribution of ^{irf)- 
In Fig. 3, the mode of the distribution is plotted as 
a function of irf/Iq- Fot irf/Io < 0.125, the bifur- 
cation amplifier is always in state 0, <p is constant and 
assigned a value of Odeg. As the drive current is in- 
creased above irf/Io = 0.125, thermal fluctuations are 
sufficiently large to cause transitions to the 1 state. In the 
region between the two dashed lines at Irf/Io = 0.125 
and irf/Io = 0.160, <f> displays a bimodal distribution 
with peaks centered at and 74 deg with the latter cor- 
responding to the amplifier in the 1 state, as we have 
demonstrated previously The dotted line in Fig. 3 
is the average reflected signal phase (<fr) . When i rf / Iq is 
increased above 0.160, the system is only found in state 1. 
In the decreasing part of the irf ramp, the system does 
not start to switch back to state until irf/Iq = 0.065. 
The critical switching currents Ir for the — > 1 transi- 
tion and I§ for the 1 — > transition, calculated from nu- 
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merical simulations to treat the inductance of wire bonds, 
are denoted with ticks in Fig. 3, and are in good agree- 
ment with experiment. The hysteresis I B < I B is a conse- 
quence of the asymmetry in the escape barrier height for 
the two states. Thus, the — ► 1 transition at i B p = Ib 
is nearly irreversible, allowing the bifurcation amplifier 
to latch and store its output during the integration time 
set by the sensitivity of the follower amplifier. 

To determine the sensitivity of the bifurcation ampli- 
fier, we have characterized in detail the switching in the 
vicinity of the — ► 1 transition. We excited the sys- 
tem with two different readout pulse protocols. In the 
first protocol, the drive current was ramped from to its 
maximum value in 40 ns and was then held constant for 
40 ns before returning to 0. Only the final 20 ns of the 
constant drive period were used to determine the oscil- 
lation phase with the first 20 ns allotted for settling of 
the phase. Histograms taken with a 10 MHz acquisition 
rate are shown in Fig. 4. In the upper panel, the two 
peaks corresponding to states and 1 can easily be re- 
solved with a small relative overlap of 10~ 2 . The width of 
each peak is consistent with the noise temperature of our 
HEMT amplifier. In this first method, the latching prop- 
erty of the system has not been exploited. In our second 
protocol for the readout pulse, we again ramp for 40 ns 
and allow a settling time of 20 ns, but we then reduce 
the drive current by 20% and measure the reflected sig- 
nal for 300 ns. In that latter period, whatever state was 
reached at the end of the initial 60 ns period is " latched" 
and time is spent just increasing the signal/noise ratio of 
the reflected phase measurement. As shown in the lower 
panel of Fig. 4, the two peaks are now fully separated, 



o 



CO 
-*— » 

o 
o 




40 80 

phase $ (deg) 



120 



FIG. 4: Histograms of the reflected signal phase 4> at Irf /Iq — 
0.145. The upper histogram contains 1.6 x 10 6 counts with a 
measurement time r m = 20 ns. The lower panel, taken with 
the latching technique, has 1.5 x 10 5 counts with a measure- 
ment time Tm = 300 ns. Data here has been taken under 
the same operating conditions as in Fig 3. The dashed line 
represents the discrimination threshold between the and 1 
state. 
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FIG. 5: Switching probability curves at T = 280 mK as a 
function of the drive current irf- The discrimination power 
d = 0.57 is the maximum difference between the two curves. 
The measurement protocol is the same as shown in the upper 
panel of Fig. 4. Though the two curves differ by approxi- 
mately 1% in 7q, they are shifted by 6% in drive current. 



with a relative overlap of 6 x 10~ 5 allowing a determi- 
nation of the state 1 probability with an accuracy better 
than 10~ 3 . This second protocol would be preferred only 
for very precise time-resolved measurements of Iq or for 
applications where a low-noise follower amplifier is im- 
practical. 

A third experiment was performed to study the state 
1 switching probability Pq—>i (iRF) for different values of 
the temperature T and Iq, the latter being varied with 
a magnetic field applied parallel to the junction plane. 
Using the first readout protocol and the discrimination 
threshold shown in Fig. 4, we obtain the switching prob- 
ablity curves shown in Fig. 5. Defining the discrimi- 
nation power d as the maximum difference between two 
switching probability curves which differ in Iq we find 
that at T = 280 mK, d = 57% for AI /I = 1% - the 
typical variation observed in a superconducting charge- 
phase qubit • The switching probability curves should 
shift according to {AI B /I B )/(AI /I ) = 3/ (4a) - 1/2 + 
0(l/(aQ) 2 ), which for our case takes the value 5.6. In 
Fig. 5, the curves are shifted by 6%, which agrees well 
with this prediction. For the case of the DC current bi- 
ased junction, similar curves would shift only by 1% since 
the switching current is Iq itself. Comparable discrimi- 
nation power using DC switching has only been achieved 
in these devices at T < 60 mK. As the temperature is 
increased, the switching probabilility curves broaden due 
to increased thermal fluctations and the disciminating 
power decreases: at T = 480 mK, d = 49%. 

We finally determined the escape rate 

IW(iHF,Jo,T) = u> a /2n exp (-MJ* y Jk B T) as a 

function of irf by measuring the time dependence of 
the switching probability, using a method previously 
applied to the determination of the static switching 
rates to the voltage state ^8|. Here oj a is the attempt 
frequency and the barrier height can be approxi- 
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ramp (40 ns) and settling period (20 ns), the reflected 
signal phase was extracted every 20 ns for a duration 
of l/is. By repeating this measurement, we generated 
switching probability histograms which we analyzed 
as Po^i (t) = 1 — cxp(— ro->it). In parallel with 
these dynamical switching measurements, we have run 
static switching measurements with the same analysis, 
yielding the escape temperature T s e / C . We have found 
that T s e t sc exceeded T by 60mK, an effect we do not 
fully understand yet but which might be attributed 
to a combination of insufficient filtering in our RF 
amplifier line outside the measurement band and to 
high-frequency non-equilibrium vortex motion in the 
superconducting films of our device. Analyzing the 
dynamical switching data with T s e / C in place of T, we 
extract a value of u° dyn = 10.7 K from the T = 280mK 
data with Iq — 1.17 /iA while the calculated value 
keeping higher order terms in 1/aQ is lt^° ° = II. OK. 

With the Josephson bifurcation amplifier operating at 
T s e t sc = 340 mK, it is possible to resolve with a sig- 
nal/noise ratio of 1 a 1% variation in Iq in a total 
time < 80 ns, corresponding to a critical current sen- 
sitivity of Sj^/Io = 2.8 x 10- 6 Hz~ 1/2 . This value is 
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in agreement with the approximate prediction S/ cx 

[diAU^/k^/dlo]- 1 ■ tU 2 - This formula shows that 
the sensitivity is independent of Q, in contrast with the 
sensitivity obtained in the linear regime [lflj . Operat- 
ing the junction near a bifurcation point has the advan- 
tage that the sensitivity is limited only by temperature 
and not by the measurement bandwidth. The advantage 
of the bifurcation amplifier over SQUIDs [20j resides in 



its extremely low back-action. Since there is no on-chip 
dissipation, the only source of back-action is the junc- 
tion phase fluctuations induced by our matched isolator 
load at T — 280 mK. Both numerical simulations and 
analytical calculations predict a spectral density S$ ~ 
10~ 12 rad 2 / Hz. In our system, the generalized back- 
action force is /(£) = — h/2e cos 6(t) and we calculate 
a noise temperature Tjv = y^SfSig/kg < 10 mK, con- 
siderably lower than that of resistively shunted SQUIDs 
where Tjv ~ T at comparable temperatures [5l| . Finally 
the bifurcation amplifier does not suffer from quasiparti- 
cle generation associated with hysteretic SQUIDS |J| and 
DC current-biased junctions |2| which switch into the 
voltage state. Long quasiparticle recombination times 
at low temperatures limit the acquisition rate of these 
devices while the recombination process itself produces 
excess noise for adjacent circuitry |22| . 

In conclusion, the Josephson bifurcation amplifier is 
competitive with the SQUID for applications where low 
noise temperature is required. Its speed, suppression of 
on-chip dissipation, and latching make it ideal for the 
readout of superconducting qubits. At temperatures such 
that T^yn — 60 mK, the discrimination power would be 
greater than 95%, hence permitting stringent tests of 
Quantum Mechanics, like the violation of Bell's inequal- 
ities. 
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